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We present a quantitative study of energy flow away from jets by numerically solving the evolution 
equation derived by Banfi, Marchesini and Smye (BMS), and apply the result to two processes at the 
LHC: Discriminating high-pt jets originating from decays of heavy electroweak bosons from the QCD 
background, and the survival probability of BFKL-initiated dijet rapidity gaps. As a byproduct, 
we find a hidden symmetry of the BMS equation which is a remnant of conformal symmetry. 
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I. INTRODUCTION 



Q . As we are entering the era of the Large Hadron Collider (LHC) , the importance of jets appears more prominent than 
' ever before. Since many of the signatures of new physics are expected to be observed via hadronic jets in disguise, 
every bit of information from jet measurements is potentially of great significance and should be carefully analyzed 
' in light of the QCD expectations. [See recent reviews and references therein.] However, efforts in this direction 

often encounter difficulties associated with energy flow — the transfer of energy or the transverse momentum away 
f~i ■ from hard jets due to the multiple emission of soft gluons. The interjet energy flow affects many observables either 
Oh' directly or indirectly, and tends to reduce the precision in their experimental measurements. For instance, a small 
uncertainty of 1% in the jet energy can lead to a 10 % uncertainty in the jet cross section at pt ^ 500 GeV 3]. 
' Normally, the soft radiation responsible for energy flow is encoded in event generators where it is no longer cleanly 
, separable from nonperturbative effects such as the underlying event and hadronization effects. Still, one would like 
' to fully understand at least the perturbative part of the radiation at a quantitative level, and indeed there has been 
substantial progress in theory along this line over the past decade. Given such progress, it would be very interesting 
^ ' if one could treat energy flow as a useful tool for discovering novel phenomena at the LHC, rather than deeming it a 
\^ I nuisance. 

. In perturbative QCD, observables related to the interjet energy flow typically involve two hard scales E ^ Eout ^ 
' Aqcd where E and Eout are the jet and interjet energy scales, respectively. In the weak coupling expansion of such 
^— ^ , observables, logarithmically enhanced terms of the form (as In E/ Eout)"' appear due to the miscancellation of real 
0^ ■ and virtual contributions. These logarithms fall into two classes which are equally important at realistic energies in 
' collider experiments: (i) The Sudakov logarithms arise from the direct emission of soft gluons from the primary hard 
partons (jets), and are thus sensitive to the antenna structure of a given process. They can be exponentiated by a 
sophisticated resummation procedure [1, 0,0- (ii) The nonglohal logarithms, first pointed out in are generated 

by soft, large-angle emissions from secondary (and ternary, etc.) gluons. These gluons multiply exponentially and 
form a cascade in the interjet region. One then has to consider the coherent emission of soft gluons from the whole 
cascade whose structure itself is determined by the previous soft radiation. Because of this complexity, the nonglobal 
' logarithms do not exponentiate, but their resummation required somewhat unconventional strategies. 

To date, two equivalent methods to resum the nonglobal logarithms have been proposed. Initially, Dasgupta and 
Salam developed a Monte Carlo simulation code to actually generate the cascade on a computer 0,1^. This approach 
was followed by several works where it was primarily used to calculate cross sections which involve rapidity gaps 
0) EH J as well as to appraise the accuracy of event generators regarding nonglobal observables [l^ . Alternatively, 
Banfi, Marchesini and Smye 13] (BMS) have reduced the problem to solving a nonlinear integro-differential equation. 
This has paved way for the remarkable correspondence between the nonglobal logarithms and the BFKL logarithms 

[a [a [ME [a. 

In this paper we adopt the second approach and perform a detailed study of the BMS equation. We then apply 
the results to two processes at the LHC where energy flow plays an important role. Firstly, we suggest the possibility 
of using energy flow as a discriminator of high~pf jets initiated by highly boosted electroweak bosons from the QCD 
background. As recently stressed in 19], such jets may result from the decay of TeV-scale new particles, so finding 
methods to identify them is an urgent task at the dawn of the LHC. Secondly, we consider the perturbative survival 
probability of large rapidity gaps created by the BFKL Pomeron exchange. We show how to include the effect of 
the finite jet cone size in the gap cross section calculated in the BFKL framework. It should be said that in these 
phenomenological applications, we do not intend to present a complete treatment of the problem. Rather, we shall 
focus on the perturbatively calculable part of the process involving energy flow which, with more work, can be 
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successfully interfaced with other details of the collision event. 

In Section II we briefly review the BMS equation and show that in a certain case it admits a 'bonus' symmetry 
which is essentially a remnant of conformal symmetry. In Section III we solve the equation numerically and present 
the result in various forms that are useful for later purposes. We then devote Sections IV and V to the applications 
mentioned above. 



II. BMS EQUATION AND ITS HIDDEN SYMMETRY 
A. The equation 

In [l^, Banfi, Marchesini and Smye (BMS) have proposed a method to quantify the efficiency of energy transfer 
away from hard jets summing all the single logarithmic terms in the large-iVc approximation. To explain this, consider 
back-to-back jets with invariant mass ^/s — 2E produced in e+e^ annihilation (Fig. [1]). Define the 'in' region Cm as 
a pair of cones with opening angle 9in pointing to the thrust axes. The complementary region is denoted as Cout- Let 
P be the probability that the total energy (or the total transverse momentum with respect to jets) radiated into Cout 
is less than Eout 

J2 E^<Eout. (1) 

It has been shown in that P obeys a nonlinear integro-differential equation in the regime E ^ Eout, provided 
one generalizes its definition to allow the two jets to point to arbitrary directions (17^, ^p) inside the two cones as in 
Fig.[i;a,b). The result is 

JCaut 47r (1 - cos 6Iq^)1 - cos 6*^/3) 



(Prij 1 — cos 6*0/3 

c-^ 4:TT (1 — cos 6*0,^) (1 — cos 0^/3) 



(^Prina,n^)Pr{fl^,ni3) - (17^^ , 17/3)) , (2) 



where the ^,;„-dependence is implicit and the solid angle integrations are restricted as indicated. The evolution 
parameter r reads 



as In , (3) 

^out 



in the fixed coupling case («<, = asNc/n), and 



r=^Un^^^^^ln(l + 2basiEout)ln^) , (4) 
TT 2b as{E) 27r6 V Eout J ^ ' 

with b = (llA^c — 2n/)/127r, in the running coupling case. The initial condition is Pt=o ~ 1 for any points (J^q, i^/s), 
and one also has that Pr{^,fl) = 1 for all values of r (no radiation from a dipole with zero size). The important 
feature of the equation ([2]) is that it resums all the single-logarithmic contributions simultaneously. The Sudakov 
logarithms are included in the first term on the right hand side and the nonglobal logarithms in the second, nonlinear 
term. 

As already noted in [13], the equation can be broadly generalized to other hard processes such as hadron-hadron 
collisions by changing the definition of Cm- Namely, the two cones need not be pointing back-to-back as is relevant to 




3 




e~^e~ annihilation in the center~of-mass frame, but their relative directions and sizes can be suitably chosen for the 
process of interest. [Though this may cause complications when solving the equation numerically.] It is even allowed 
that Cin consists of more than two cones, or only of one cone as in Fig.[2l[c). As a matter of fact, we have found that 
this single-cone configuration is a particularly interesting case both from mathematical and phenomenological points 
of view, and therefore it will be our main focus in the following. 



B. Hidden symmetry 

Unfortunately, the nonlinear equation ^ is too complicated to be solved analytically, though a universal feature 
("geometric scaling") arises in the large-r region for which some analytical insights can be given flS*!. In realistic 
collider experiments, there is a rather severe restriction t < 1 ~ 1.2 , and in this regime the equation has to be studied 
numerically. Still, here we show that in the single-cone case there exists a hidden symmetry which puts a strong 
constraint on the solution. For this purpose, it is convenient to employ the exact correspondence p7lll8|| between the 
interjet soft gluon cascade and the BFKL dynamics [20, 21]. The equation ([2]) defined on a two-sphere with the 
coordinates O — {9,(j>) can be mapped onto an equation on a two-dimensional transverse plane x = (x^,x^) via the 
stereographic projection (see. Fig. ^ 

9 6 
= tan — cos (f> , x^ — tan — sin ^ . (5) 

In the single-cone case, we take the origin of the polar coordinate system in the direction of the cone axis and obtain 

drPriXa^Xp) ^ -fapPriXa^Xis) + f ^T"^ 7^ ^^-^^2^^^ (-Pr (S'a , ^7 )^V (^7 , ^/j) - (^a , )) , (6) 

J\x.y\<ri„ l-^" Xj) (Xj Xp) \ / 

where 

cPx^ {xa — xpY 

'\x-,\>rin 27r (Xa ~ Xj)'^{x^ — Xfi) 



The solid angle integration inside/outside the cone has been mapped onto an integration inside/outside a disk with 
radius 

Tin = tan ^ . (8) 

Let us rescale x/r^^ — > x so that ri„ = 1. If one relaxes the restriction \x^\ < 1, © is identical to the Balitsky- 
Kovchegov (BK) equation for the dipole S-matrix SrixaiXp) which describes the gluon saturation in high energy 
scattering [22l. 1231] . As is well-known in that context, the integration kernel 

2tT {Xa — Xj)'^{Xj — XpY ' 

is invariant under conformal transformation which forms the group SL(2,C) 

z = x^+ix' ^z' ^ ^1±1 , aS-l3-f^l. (10) 
7z + 
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Introduction of the boundary ri„ = 1 breaks conformal symmetry down to the subgroup SU(1, 1) ~SL(2,M). This 
maps the interior of the disk onto itself, and is defined by the foUowing transformation 

z^z'^f^. \a\^-W^l. (11) 

pz + a 

In the context of the BK equation, conformal symmetry is of limited use because it is broken by the initial condition. 
However, the initial condition of the BMS equation does have the SU(1,1) symmetry trivially because PT-=o(^^a, ^p) = 
1 for all ria,/3- A disk which has the isometry group SU(1,1) is known as the Poincare disk. The invariant measure of 
the distance (the chordal distance) between two points {Xa,xp) on the Poincare disk is (see Appendix A) 

d'{x^,x,)^ (12) 

(1 - xijii-xjj) 

Since both the equation and the initial condition are SU(l,l)-invariant, the solution Prixa, x/^) must also have this 
symmetry which means that it is a function only of (P(xa, xp) (and t), or equivalently, only of the geodesic distance 



l{Xa,Xi3) = cosh ^(1 + 2d^(xQ,X^)) = 2cOsh ^ y^l + (P{Xc,,Xf3) . (13) 

In order to demonstrate the usefulness of this observation, consider the integral This can be easily evaluated 
if one sets xp = 0. The result is, restoring rm, 

/ao = ^ In ^— . (14) 

2 1 — ^ 



xl 



Using the relation 

d'(a?a,0) = (15) 
one can write 

/aO = ^ln(l + d'(fa,0)). (16) 
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Since /^^ is a function only of <P{xa,Xj3), the result with generic 7^ is simply given by 
/„^ = lln(l + d2(f„,x^)) = ilnfl ' ^-(^"-^/')' ^ 



2 2 (r?„ - xiM^ - xj) J 

I sin^6l,„(l-cos6l„;3) \ ^ 

2 \ 2(cos — cos 6'i„)(cos 6'/3 — cos ^in) / ' 

where in the second line we switched back to the original sphere problem using the stereographic projection. It 
requires a considerable amount of work if one tries to get the same result for fa/3 by directly evaluating the integral 
d?]) with X0 ^ 0. 

Similarly, the full solution PT{xa, xp) can be obtained from Pr{xa, 0) which is a function only of \xa\/rin 



PAXarO)^Pr[^^-Pr . (18) 



Prixa, Xfj) for generic xp ^ is. given by 



Pr{Xa,X,) = P. LM^^fe^ - Pr i ) . (19) 



1 + (f{Xa,Xi3) I \\1 - ZaZpl 

In practice, when solving ([6]) (or ([2])) numerically, one registers the values of Pr for all (discretized) points (xa^xp) 
at each step of iteration. This can be done straightforwardly without really caring about the constraint p^ . The 
advantage of is that a single plot of the function Pt- (xq, , 0) thus obtained tells the value of Pr for an arbitrary 
point {xa, Xfj). 

Going back to the sphere problem, psp implies that 



/ tan — \ 

Pr{n, 0)^Pr{ = P.(e''--") , (20) 

\^tan^y 

where the pseudorapidity variable 77 is defined as usual 

?7==lncot-. (21) 

(1^0)1 shows that, if one of the jets is along the cone axis, P becomes a function only of the relative rapidity between 
the other jet and the cone edge. This is of course a consequence of boost invariance. In the x coordinate, it has a 
very simple meaning as the dilatation symmetry x — > ex. Once the cone size is fixed, the boost invariance is lost, 
and naively one would expect that the only symmetry of the solution Pr{^lon^[3) would be the trivial 0-rotation. 
Remarkably, however, a part of conformal symmetry of the soft emission probability survives, and this reduces the 
four degrees of freedom {^lajflp) to one. Note that the SU(1,1) symmetry no longer exists if there are two cones 
forming Ci„. However, we shall see that the difference in the numerical solutions with and without the backward cone 
turns out to be quite small. 



III. NUMERICAL RESULTS 

In this section we present the numerical solution of the equation ^ extending the initial result in . In doing 
so, it is useful to factor out the Sudakov contribution 

Prina,np) = e-^f''^grina,np), (22) 

where fa/3 is given by (fT7|) in the single-cone case (Fig. Etc)), and by Eq. (C.l) of in the case where two cones 
are pointing back-to-back (Fig. [2)Ja,b)). The effect of the non-global logarithms is then contained in gr- At small 
T <C 1, only the Sudakov contribution is important, whereas at large r 3> 1 the nonglobal contribution dominates in 
the sense that e""^-^ ^ 3t- In the phenomenologically relevant range r ~ 1, the two contributions are comparable. 
Note that in the single-cone case, / and gr are separately SU(l,l)-symmctric. 
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FIG. 4: Solid curves: the numerical solution of Eq. ((2| in the two-cone case. Dotted curves: the Sudakov contribution. Four 
curves correspond different values of r: From top to bottom, r = 0.6, 0.8.1.0, 1.2. 



Let us first consider the two-cone case. We fix Op — 0, Oin = 7r/3 and plot Prida = 0, 0) for four values of r between 
0.6 and 1.2 in Fig. [H The dotted curves denote the Sudakov contribution alone. The result shows the expected 
behavior; for small 9 corresponding to a dipole (dijet) with small opening angle, energy flow is suppressed due to the 
QCD coherence and P is close to 1. On the other hand, P is significantly less than 1 for dijets forming a large angle 
6 > 2tt/3. One also sees that the nonglobal contribution becomes more important as r gets larger. 

Next, we consider the single-cone case and show in Fig.[n]the 'master function' (pOI in the same range 0.6 < r < 1.2 



as a function of 
this case reads 



tan 9/2 
tan 9i„/2 



(left) and of 77 — ryi„ (right). Again the dotted curves represent the Sudakov factor which in 



-Tf _ 



= 1- 



tan 



tan 



2 9i. 



(23) 



In fact, numerically the result in Fig. [5] is very close to the left branch of Fig. 21 that is, the difference between 
Figs. [UJb) and[2Kc) is tiny, less than 1% even for 9in as large as tt/S. This means that when the two jets are in the 
same cone the effect of the backward cone is negligible, and practically one has the SU(1,1) symmetry to a very good 
approximation. 

Finally, we consider the probability that the interjet energy fraction 



E 



(24) 



is exactly Cout = Eout/E. This may be obtained by differentiating P^- with respect to Eout- For simplicity, we work 
in the fixed coupling case and denote the distribution of Cout as W{eout)- By definition. 



W{e) de = Pr 



or equivalently. 



W{eout) = 



-Pr 



-ase 



r,^,dPr_ 



(25) 



(26) 



The function W{eout) with = 0.17 is plotted in Fig. [H] for the back-to-back case (left) and the single-cone case 
(right); the latter is obtained from (|20p with a particular value tan |/ tan ^ = 0.46 to be relevant later. W is 
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FIG. 5: Numerical solution of Eq. ([2]) in the single-cone case. 



sharply peaked in the small-Cout region {eout ^ ^s) where the soft approximation is rcUable, and the suppression in 
the 'large'-Cout region (cout ~ Q^s) is stronger in the single cone case as expected from the coherence effect. Though 
not shown in the figure, W actually starts to decrease as Cout becomes extremely small {eout < 10~^). This can be 
understood from the asymptotic large-r analysis in [1, [l^ . 
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FIG. 6: The probability distribution of interjet energy fraction W{eout)- Left: Two-cone case with back-to-back jets. Right: 
Single-cone case. 

Fig. [6] shows that in the majority of events the interjet energy flow is very small Eout ^ ctgE. These typical events 
do not contribute to the average of Eout which is dominated by rare events with Eout ^ oisE. Indeed, in the present 
context it is tempting to define the RverRge of Equi cls 

{eout) = / deout W{eout) eout = - / dr e"* — = 1 - — / dr e-^/"=P, , (27) 
Jo Jo cIt a, Jo 

where in the last equality we have integrated by parts. Clearly, the integral is dominated by the region t < as 
corresponding to > Eout ^ cXsE. The soft approximation is not reliable for such large values of Eout, one has to 
use instead the full splitting function to compute {Eout) (see, e.g., [24|). Thus, although the probability distribution 
W{eout) is meaningful at small eout <C 1, it is not entirely legitimate to compute the average of eout using (|27|1 . 
Nevertheless, the calculation is simple enough and can be done analytically since it is dominated by the Sudakov 
contribution. Moreover, the result may be useful to infer the qualitative feature of the differential spectrum of Eout 
in the interjet region. We outline such calculations in Appendix B. 

Having discussed the detailed properties of the probability distribution P^, we now turn to its phenomenological 
applications in collider experiments. 
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FIG. 7: Left: Massive weak boson jet at midrapidity. Middle and Right: Massive quark jet together with an accompanying 
jet ('color neutralizer'). 



At the LHC, highly-boosted heavy electroweak bosons {W, Z, Higgs) can serve as a signal of new physics. For 
example, they arise as decay products of TeV-scale new particles in certain extensions of the standard model (see, [Toj 
and references therein) . In order to maximize the potential of discovery, clearly it is desirable to be able to identify 
these events in hadronic (as well as leptonic) decay channels. Then the issue arises as to how one can experimentally 
distinguish, preferably event-by-event, massive jets originating from boosted weak bosons (signal) from QCD jets in 
the same mass range (background) initiated by light quarks and gluons [1^, [2^ . Ref . [1^ looked into the substructure 
of jets for this purpose and showed that it has different characteristic behaviors depending on the progenitor. [For 
related works, see [23, HI, HI, HO, HII-] Somewhat complementary to this, here we investigate the pattern of energy 
flow outside the jet cone as a possible discriminator of weak boson/QCD jets. Our goal in this section is not to make a 
practical proposal that is readily useful in experiments, but rather to give the first quantitative study of the difference 
in energy flow between the two types of jets which will lay the foundation for future work. 

Fig.[7]shows a 'two-pronged' jet at midrapidity originating from a high-p^ boson (say, Z^) of mass M pt together 
with those from a quark with virtuality M . The opening angle between the two primary decay particles inside the 
jet is bounded from below 



in both cases, with the lower bound Oap ~ 2M/pt being the most probable configuration. Away from this peak, the 
two types of jets have different particle distributions in [2^ . 

In the case of a weak boson jet, radiation from the (color-singlet) qq pair is strongly suppressed due to the QCD 
coherence. Let us quantify this statement by computing the probability P-r discussed in the previous section. We 
assume that the jet is produced around midrapidity « and the qq pair is symmetrically emitted so that fJo, — [9, 0), 
Q,p = {6,tt) with tand = M/pt, where the polar angles are measured with respect to the triggered jet axis. We then 
identify the jet cone with Ci„ (i.e., the single-cone case) by taking Om ~ R where R = y/ l^rf + Ai/)^ is the usual jet 
cone radius. [Around midrapidity, A9 k, Arj for small R = Arj.] Using the SU(1,1) symmetry p9)) . one can show that 
the configuration fi^ = (^^,0), il/j = {0,tt) is equivalent to the configuration 



IV. HIGH-pt JETS FROM HEAVY ELECTROWEAK BOSONS 




(28) 



(29) 



where 0' satisfies 
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(30) 
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This means that 
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FIG. 8: Red curve: The probability as a function of -^^jjr- Green curve: Sudakov contribution ((33 



namely, the probability distribution can be obtained from the master function Pt{x) merely by the change of variables. 
We take E = pt = I TeV and Eout = 10 GeV with as{Eout) ~ 0.18 in jH) to get r w 0.6. In Fig. [3 we have plotted 
Pr=o.6 as a function of 



tan^ ^ ptR 
tan I ^ ~M" 



(32) 



The upper curve is the Sudakov contribution which in this case reads (cf., (|23p and ([50)) ) 

(33) 



tan ■ 
^- : 



^ + tan2 f / / \ tan2 !^ ~^ ^ , 

Let us take, as an illustration, M — 100 GeV so that ptR/M = IQR. Then Fig. [5] shows that, if we take R — 0.4, 
with 85 % probability energy emitted outside the jet cone is less than 10 GeV. This is only 1% of the jet pt and 
is virtually indistinguishable from the underlying event and pileup contributions (assuming that the typical energy 
scale of the latter is around this value at the LHC). The situation is similar for jets with higher pt values. Even with 
Pt = 4 ~ 5 TeV, T can reach only up to 0.7 and the function (PT|) is somewhat decreased overall. However, this is 
more than compensated by the effect due to the decrease of the qq pair opening angle (if M is kept fixed) , and again 
one finds that with very high probability > 90%, i^otjt is less than 10 GeV. On the other hand, one can increase t 
appreciably if Eout is decreased, though 10 GeV might already be a bit too small value to choose in high-luminosity 
measurements at the LHC. The probability distribution of Cout — Eout/pt is shown on the right hand side of Fig. [5] 
in the previous section, albeit in the fixed coupling case. As a matter of fact, this figure has been obtained for the 
particular qq configuration above {pt/M — 10, R = 0.4). 

Next we turn to the quark jet. Clearly, energy flow is more efficient in this case. Due again to the QCD coherence, 
radiation from the qg system with 9af3 <C 1 is effectively that from the parent quark. If this quark is created via 
the qq — > qq hard scattering with one-gluon exchange, the compensating color is carried by the outgoing jet in the 
backward direction (Fig. [71 middle). In the qq qq case it is carried by one of the incoming partons (Fig. [T] right). ^ 
We thus consider the two-cone configuration with 0in = R and take Oap = tt and 7r/2 in the two cases, respectively. 
However, solving the BMS equation for the latter case is technically difficult because it involves complicated angular 
integrations. Rather than doing this, here we give a simple estimate of Pt by noting that, to lowest order in the 
small parameter R, the configuration in Fig. [T) (right) is obtained by boosting the back-to-back configuration with 



There are subleading contributions in Nc where the color flow assignment is interchanged between the two cases. 
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9in = V^R in a direction orthogonal to the jet axes with velocity v = l/\/2. Of course this artificially modifies certain 
components of the particles' four-momenta up to a factor 7 — l/Vl — = v^, resulting in a small change in t 
mainly via the running of the coupling. [The ratio pt/ Eout is less affected by the boost.] But in the present study we 
ignore this change. 

The results are plotted in Fig. [9] as a function of R. In both cases, Pr is significantly smaller, by a factor of about 
five, than in the weak-jet case. With very high probability (80 ~ 90% when R = 0.4) the energy radiated outside the 
jet cone is greater than Eout = 10 GeV, though it is typically smaller than the average value {Eout) ^ ctsPt ^ 0(100) 
GeV as can be inferred from the left hand side of Fig. [6l 

We have thus demonstrated the striking contrast in the amount of energy flow between jets initiated by colorless 
and colored particles. The five-fold difference in the probability Pr at first looks a powerful criterion for the event- 
by-event identification of jet types. [Note that the average energy is dominated by rare events, so it is not suited for 
the event-by-event identification.] However, in realistic experiments several other effects are expected to reduce this 
difference. Most importantly, in the above we have dealt with an idealized situation where there is only one dipole 
containing the triggered jet, and said nothing about the presence of other dipoles in the process some of which are 
shown in Fig. [7] in dotted lines. The radiation from these extra dipoles is presumably not much smaller than that 
from the triggered quark jet, and clearly one has to come up with a method to minimize these backgrounds. A step 
in the right direction would be to look at the angular region just outside the triggered jet cone where actually most 
of Eout is concentrated.^ This is briefly explained in Appendix B and will be studied further elsewhere. 




FIG. 9: Radiation outside the QCD jet as a function of the cone radius. The left (right) hand side corresponds to the middle 
(right) figure in Fig. [T] 



V. BFKL DUET CROSS SECTION WITH PERTURBATIVE GAP SURVIVAL PROBABILITY 

Our second example is the production of dijets separated by a large rapidity gap A77. Such 'jet-gap-jet' events have 
been reported at the Tevatron some time ago [32, . ,33] and are certainly an interesting channel to study at the LHC 
given the larger pt values attainable by forward jets. Two approaches to study jet-gap-jet events in perturbative 
QCD are commonly used: 

(i) In the BFKL approach [si, [H, [H, [13, El, El] , the exchange of the BFKL Pomeron [13, HH in the i-channel 
naturally generates a rapidity gap, which should be the dominant process in the limit of a large gap Ayy 3> 1. At 
the Tevatron, gap events are identified by requiring no activity (above the threshold) in the central rapidity region 
\ri\ < 1. Yet it seems more natural to define a 'gap' in the true sense of the word, to be the region between the edges 
of jet cones so that the actual gap size is A77' — Arj — 2R where Arj is the rapidity difference of dijets. This is desirable 
also from a theoretical point of view because in this way one can efficiently suppress the contribution from color-octet 



^ Related to this point, we note that the dipole number (not energy) distribution emitted from a boosted qq pair can be computed exactly 
including the nonglobal logarithms Il8|| . 
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FIG. 10: Radiation from two independent color dipoles. 



exchanges so that the signal of the BFKL Pomeron becomes more visible. However, jet cones filled with soft radiation 
are simply absent at the level of the BFKL cross section, though an event generator may be used to fix this. 

(ii) In the factorization approach 0, [^, [l^, [4^ , one defines a gap event in terms of energy flow and starts out 
with the one-gluon exchange. One then dresses the amplitude with soft gluons via the renormalization group with 
due respect to the constraint in energy flow into the gap region. The BFKL Pomeron does not arise automatically in 
this approach, rather, it has to be added by hand avoiding double counting [43|. The nonglobal logs are not included, 
but their effects may be estimated separately ^10]. 

Here we employ the BFKL approach and discuss how to include the effect of soft radiations and of the finite jet cone 
size. We define gap events in terms of energy flow (as in the factorization approach) and require that the total amount 
of energy between the edges of jet cones is less than Eout- Events with a perfect gap would correspond to Eout = 0, 
but in order to use perturbation theory one should require Eout ^ ^qcd- We take, somewhat optimistically, Eout — 1 
GeV and 2 GeV having in mind low-luminosity measurements.'^ 

Specifically, consider pp or pp scattering with center-of-mass energy s = (2i?)^. We parameterize the momenta of 
the underlying partonic process ab — > 12 as (see. Fig. (TUl) 

Pa = {XaE,0,0,XaE), (34) 

Pb = [xbE, 0,0, -XbE), (35) 
pi = {pt cosh T]i, pt, Pt sinhr]i) , (36) 
P2 = (ptCOsh 772,-pt,Ptsinh772) . (37) 



The relative rapidity of the two jets is A77 ~ rji — rj2, while the average is 77 
The cross section in the BFKL approximation is given by {t « —Pt) [13] 



_ vi+m 
2 



da da r , r , ^ da?'?-'?'? , , 

^a^b-] 1 T7 = XaXbJeff{Xa,t)Jeff{Xb,t) — , (38) 



dAT]dfidt dXadXbdt ' ' ' ' dt 

where the effective parton distribution is {Cp — {N'^ -~ 1)/2N, 



f,ff{x,t) ^ q{x,t) + q{x,t) + -4-g{x,t) , (39) 



and at the leading logarithmic level, 



> (40) 



dt ttN^ t^ \J (i/2 ^ 1/4) 

with oj{i') = 2Re [^/'(l) ~ tp (^^ + ii^)]- Attempts to include the higher order conformal spins and the next-to-leading 
logarithmic effects can be found in HI) H^l • 



^ Note that, as already implicit in the previous section, in our definition Eout includes only the perturbative radiation from the hard 
scatterers, that is, the contributions from the underlying event and pileups are subtracted. These are partly taken into account by the 
nonperturbative survival probability S introduced later. 
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We now implement the survival probability of the gap against soft radiations. It is well-known that the non-forward 
BFKL amplitude coupled to quarks contains the Sudakov factor 

exp(^-a,Ar;lng^ , (41) 

where kt <C pt represents the loop momentum in the soft region, arising from the exchange of Reggeized gluons in the 
t-channel. This limits emissions into the gap region, but apparently, the amplitude is sensitive to the infrared region. 
However, as explained in [4lj such sensitivity disappears in the Mueller-Tang formula [s^. This renders the partonic 
cross section (PH)) truly dominated by short distance physics at the scale ^/\t\ = pt, thereby allowing one to treat it 
nearly on the same footing as, say, the one-gluon exchange contribution. Once this has been done, however, there 
arises additional possibility to fill the gap by soft radiation from color dipoles (paPi) and {pbP2)- [Since the BFKL 
exchange is color-singlet, color flows as a — > 1 and h — s- 2.] This is not included in the BFKL approximation since the 
corresponding diagrams are not enhanced by powers of A77. Still, they are enhanced by powers of hipt/Eout and can 
be taken into account by the survival probability Pr- 

In the large Nc approximation, radiations from the two dipoles (paPi) and {pbP2) are independent. We thus modify 
,s 

da ^ ^ (^Q-qg-'QQ 

— — — = XaXbfeff{Xa,t)feff{xb,t) S Pr{^l , fta)PT{^2 , ^b) , (42) 

diXrjdrjdt dt 

where 



iV2 

feff{Xa,t) = q{Xa,t) +q{Xa,t) + g{Xa, t)Prifll , i^a) 



F 

CI 



feffixb,t) = q{xb,t) + q{xb,t) + ■^g{xb,t)Pr{ft2,^b) ■ (43) 



S is the nonperturbative survival probability which concerns the soft interaction among the spectators (the underlying 
events) . It depends only on and has been estimated by several groups [H, IHI . The presence of an additional 
factor of Pr in front of the gluon density in (|43p is because at large a gluon is represented by a double line and 
radiates like the square of a quark jet. 
From (|^n|) . one has 

Pr{ni,ila) ^ Prie"^^^-'") = Prie-'') , (44) 

where, following our definition of a gap, we choose 77™ such that 771 — 77™ always equal to R, see Fig. [TOl Similarly, 
\V2 — vT\ — ^- 611) nieans that, once the cone radius is fixed, the perturbative survival probability depends only on 
Pt (or equivalently, r) of the jets, and not on the rapidity gap Arj in the present approximation. 

On the left hand side of Fig. [TTl we plot as a function of R for fixed r — 1.2, corresponding to pt ^ 300 GeV 
and Eout = 1 GeV in the running coupling case with as{lGeV) « 0.5. On the right hand side of Fig. [Til we fix 
R = 0.4 and plot (jU) as a function of pt using the one-loop relation between r and pt. The case with Eo^t = 2 GeV 
is shown in Fig. [121 We see that the results depend rather sensitively on the value of Eout in the few-GeV region. 
When Eout = 1 GeV the nonglobal logarithms are more important than the Sudakov logarithms. 

The dependence of Pr on jet pt is weak for pt > 50 GeV and is typically Pr ~ 0.5 when Eout = 1 GeV. This means 
that the BFKL cross section gets smaller at least by a factor P^ ~ 0.25. The gluon-gluon contribution is suppressed 
by an even smaller factor P^. On the other hand, by requiring the gap size to grow linearly with the dijet rapidity 
gap A77' = A?7 — 2R and choosing a small Eout , one can enormously suppress the color-octet contribution. This can 
be readily inferred from the right branch of Fig. |4l One sees that the typical suppression factor for the color-octet 
contribution is ~ 0.01 or less when r = 1.2. Thus the suppression factor ~ 0.25 in the color-singlet channel 
appears to be a small price to pay for a clean test of the BFKL effects. In practice, in order to get the full cross 
section the above results should be combined with an event generator paying attention to the double counting of the 
sort discussed in 1121. 



VI. CONCLUSIONS 



In this paper we have presented a detailed numerical study of the BMS equation which describes the dynamics of 
interjet energy flow at the single-logarithmic level in the large- A'^c approximation. In hard scattering processes, energy 
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FIG. 11: Left: The survival probability Pr (red line) at r = 1.2 as a function of the jet cone radius 7?. The blue line is the 
Sudakov factor e~'^-^ (cf. (|22p ) and the green line is the nonglobal contribution g. Right: The survival probability as a function 
of pt for fixed R = 0.4. 




flow is often thought of as an annoying feature which, like the underlying event and the hadronization effect, blurs 
otherwise clean predictions of perturbative QCD or signatures of new physics jsj. Contrary to this, our main motivation 
has been to demonstrate the usefulness of energy flow as a primary diagnostic tool to discover interesting physics 
at collider experiments. We have discussed two examples for this purpose: Identifying high-pt jets of hadronically 
decaying electroweak bosons among the QCD background, and suppressing the color-octet contribution relative to 
the BFKL Pomeron contribution in dijet rapidity gap events. Both processes are characterized by the large ratio 
Pt/Eout, and we have shown that the consideration of energy flow in terms of the probability Pt can significantly 
enhance the signal-to-background ratio. We also suggest that the probability Pr is better suited to the event-by- 
event identification of jets than the average energy (Eout), since the latter does not reflect the property of the majority 
of events. In future it remains to integrate the results into the full collision environment to see quantitatively how the 
nonperturbative effects affect the perturbatively calculated probability P^. 
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Note added. While this paper was being finalized, a preprint [44| appeared which also studied energy flow as a 
tool to identify new particles in a difi^erent context. 

APPENDIX A: POINCARE DISK 

The group SU(1,1)=S0(1,2) is the isometry group of the hyperbohc space H2 defined by 

Xl - Xl -Xl^l. (Al) 

Via the change of the coordinates 

^0+^2 = -, x = yX^, (A2) 

y 

the space is mapped onto the upper half plane with the squared length element 

ds^ = -dXl + dXl + dXi = ^ , (A3) 

y2 



where uj = x + iy and y > 0- This can be further mapped onto a unit disk via the transformation 

2 Mzdz 



ioj + 1 ,9 Adzdz 

(7S — 



A disk endowed with this metric is known as the Poincare disk. The invariant distance between two points (the 
chordal distance) on the Poincare disk can be immediately derived from that in the H2 space 

Using this, one can rewrite the kernel ([HI) in a manifestly SU(l,l)-invariant form 

J__dzjdzj d'^{za,Z0) 

2Tril-\z^\^rd^iz^,z^)d^z^,z^)- ^ ^> 

APPENDIX B: THE AVERAGE ENERGY 

In this Appendix, we estimate the rapidity (angular) distribution of Eout in the interjet region using (j27p . As 
already noted, strictly speaking this goes beyond the validity range of the soft approximation. But nevertheless we 
wish to record the result here for possible future applications. 

Let us begin with a dipole formed by the incoming partons as shown by dotted lines in the left and middle figures 
of Fig. [71 Since the r-integral in ((27)) is dominated by the small-r region, one may approximate Pr by its^udakov 
part. Actually, one can go one step further to include the first correction from the nong lobal logarithm [8|,1]J|. This 
gives 

where 

„ 01 o J, Li2tan4%^ ^2 Lise-^"- 

i?o = 21ncot— = 27?„, ^i = Y^ 2 12 2 ' ^ ^ 

[Li2 is the dilogarithm function.] Therefore, 

^1-^ rd.e-a-o).-«.^^l_^,%i^Erfcf , (B3) 



Pt as Jo 2asy/R{ \ 2y/R{ 
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FIG. 13: Left: The rapidity distribution of Eout emitted from the dipole formed by incoming partons. The region < pt/M = 
0.1 is excluded. Right: The rapidity distribution of Eout near the quark jet. 



where Erfc is the complementary error function. Let us then introduce a differential distribution 



{Eout)= I ^ dr^{^ 



= 2 



drj 



dEou 
dr/ 



(B4) 



such that 



dEo 



drj 



2 



dm, 



(B5) 



(|B5P is plotted on the left hand side of Fig. [13] with or without the Ri term. The distribution is flat to lowest order in 
tts, but after the exponentiation it slowly varies with rj. Note that Eout in this case should be understood as the total 
transverse momentum. To get the energy distribution one has to multiply the result by the factor 1/ sinO — cosh 77. 

Next we consider the back-to-back dipole including the outgoing quark jet, Fig. [7] (middle). This is basically the 
same configuration as above, but now we have to set the rapidity of the jet to be zero (rather than infinity). Then it 
is convenient to begin with the angular variable and approximate 9in ~ R for small R. One finds, omitting the i?i 
term, 



{Eo 



Pt 



TT-R 



de„ 



dO 



d6 = 



2q!s In cot 



1 



2as In cot 



2a,ln-+0(a2) 



(B6) 



If one takes only the first term in the a^-expansion, there is a logarithmic divergence as i? — > 0. This was discussed 
in (2^ but the coefficient here is different because we have included only the term in the splitting function which is 
singular in the soft limit. On the other hand, the all-order expression in (|B6[) is finite as i? ^ 0, and may be regarded 
as a partial resummation of the logarithms {a^ Ini?)". The resulting rapidity distribution {dcout/ drf) « {deout/d9) is 
plotted on the right hand side of Fig. [T31 Compared with the left hand side of Fig. [T31 one sees that the shape is 
similar, but the magnitude is bigger by a factor of two in the same rj range around midrapidity. 

Finally, Fig. [Til shows the rapidity (angular) distribution of Eout emitted from the weak boson jet with pt/M — 10 
calculated from (|33p . The distribution is very strongly peaked at small rj in striking contrast with the previous two 
cases. This is another way of seeing the difficulty to lose energy for a weak boson jet. See, also, [26| . 
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